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Abstract 

We calculate the pretzelosity distribution (hj-j,), which is one of the eight leading twist transverse 
momentum dependent parton distributions (TMDs), in the light-cone formalism. We find that this 
quantity has a simple relation with the quark orbital angular momentum distribution, thus it 
may provide a new possibility to access the quark orbital angular momentum inside the nucleon. 
The pretzelosity distribution can manifest itself through the sm{34>h — (ps) asymmetry in semi- 
inclusive deep inelastic scattering process. We calculate the sm{3(l)h — (ps) asymmetry at HERMES, 
COMPASS, and JLab kinematics and present our prediction on different targets including the 
proton, deuteron, and neutron targets. Inclusion of transverse momentum cut in data analysis 
could significantly enhance the sin(3(/>/i — (ps) asymmetry for future measurements. 
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I. INTRODUCTION 



The nucleoli is a composite system, so its spin not only comes from the intrinsic spins 
of the constituents, but also from the orbital angular momenta due to the relative motion 
of quarks and gluons. One of the main tasks of hadron physics is to quantitatively know 
the contribution from each component of the above items. However, on the theoretical side, 
there are disputes on how to define these quantities, and decomposing the contribution from 
spin and orbital angular momentum parts in experiments is also not clear yet. 

Conventionally, it seems easy to start from the QCD Lagrangian and follow the Mother's 
theorem, and to write the conserved angular momentum as 

^QCD = J d^xtp^^tp + J d^xip^x X (-iV)V' 
+ j d^xE" xA^ + j d^xE^'x X VA"' 

= Sg + Lg+Sg + Lg, (l) 

in which the four terms are identified as the quark spin, quark orbital angular momentum, 
gluon spin, and gluon orbital angular momentum, respectively. This definition has been 
used in many references 

HQ 

, and it is naturally considered as the generator of the space 
rotation. In Ref. 2|], the matrix elements of this defined orbital angular momentum operator 
in the light-cone representation were calculated, and the results satisfy the conservation of 
the z projection of angular momentum: 

n n— 1 

J' = Y.'' + Y.^'r (2) 

i=i j=i 

But a problem is that except for the quark spin term, the other three terms are gauge 
dependent, thus they have obscure physical meanings in common situations. 
By considering this, Ji suggested a gauge-invariant expression j3|: 

Jqcd = J d^x^^^i! + j d^x^'^x X {-iD)tp + j d^xx x {E x B) 

^ Sg + Lg + Jg, (S) 

where D = 'V — igA is the covariant derivative. As pointed in Ref. this can be experi- 
mentally accessed in the deep virtual Compton scattering (DVCS) process, where the total 
angular momentum Jq = Sq+Lg can be measured, and HERMES Collaboration has reported 



their experimental progresses in recent years [J]. But in this definition, unhke the quark 
case, there exists no gauge-invariant decomposition of spin and orbital contributions for the 
gluon angular momentum Jg. However, de Teramond and Brodsky showed recently that 
under the framework of light-front QCD, there is a separation of the dynamics of quark and 
gluon binding from the kinematics of constituent spin and internal orbital angular momen- 
tum. More seriously, the assumed quark orbital angular momentum operator does not obey 
the angular momentum algebra J x J = ij, and there is an interaction term involving also 
gluons. So it is a question whether such defined and measured quantity can be viewed as 
the quark orbital angular momentum. 

Trying to reconcile the gauge invariance and the angular momentum algebra, Chen et 
al. proposed a new form of spin sum [6| by adding a surface term J d?xV ■ [E^-^A^^^^ x f)], 
which vanishes after integration over the space, 

Jqcd = j d^xi)^'^^ + j d^xip'^xx (-iZ3pure)^ 
+ j d^xE"" X I^hys + j (i^xE^'x X VAj^i^y, 

= Sg + Lg + Sg+ Lg, (4) 

where -Dpure = V — igAp^j.^, and Apure is the pure gauge component of A, satisfying -Dpure x 
^pure = 0- A can be decomposed as A = Apure + ^phys- In this definition, all terms are 
gauge invariant and obey the angular momentum algebra. But there is a practical problem 
of lacking a clear way to measure orbital angular momentum defined so far. 

If we disregard the disputes on the definition, we can define various quantities under each 
frame of definition and investigate their properties. In Ref . j?! , the first definition mentioned 
above was used and a quantity that represents the quark orbital angular momentum in the 
light-cone formalism was obtained. It was shown that the quark orbital angular momentum 
equals the difference between the helicity and transversity, but it was not clear how to 
measure this quantity directly at that time. In Ref. j^, Avakian et al. used the bag model 
and the spectator model to deduce the result that the difference of helicity and transversity is 
(a transverse moment of) the so-called pretzelosity, a new quantity which is one of the leading 
twist parton distributions. Later in Ref. p], Pasquini et al. confirmed this statement but 
using the light-cone constituent quark model. Conventionally, the difference of helicity and 
transversity is viewed as the relativistic effects in the nucleon, so in Refs. j^, 9], pretzelosity 



3 



is considered as the measurement of the relativistic effects in the nucleon. In this paper, 
we calculate the pretzelosity distribution in the light-cone SU(6) quark-diquark model and 
reconfirm the conclusion as that in Refs. . So combining with the result already obtained 
in Ref. 7[, we suggest that pretzelosity might open a new way to access the quark orbital 
angular momentum. Of course we should remember that our definition and conclusion are 
based on the light-cone gauge in the infinite momentum frame (IMF), or equivalently, = 



gauge in the light-cone formalism 



10 



ll|, which is the conventional language for describing 



parton distributions inside the nucleon. The advantages of using light-cone formalism to 
describe the quark orbital angular moment have been also stressed in Refs. {2, 5]. 



II. TMDS AND PRETZELOSITY 

Pretzelosity, usually denoted as h^rp, is one of the eight leading twist transverse dependent 
parton distributions (TMDs). Up to the leading twist, the decomposition of the quark-quark 



correlator reads 



12 



13| 



Using the abbreviation = Tr($r)/2, we can get the traces of the correlator jl2l. \\\ 

= (6) 



N 



$'^"^^1 = + ^^^?1T, (7) 

. s. + s„ + s'J^^ki, - (8) 

Working in the light-cone gauge, we have HQ 

The pretzelosity distribution can be worked out by 

^/^^t(^,p1) = / '^^^e^^^^^^--^--^-\psy\m^<^'^i^m (10) 

where \PSy) denotes the hadronic state with a polarization in y direction. 
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Pretzelosity is chirally odd, so it must couple with a chirally odd partner in the semi- 
inclusive deep inelastic scattering (SIDIS) process (see Sec. V). It is usually considered as a 
measure of the relativistic effects in the nucleon, and a detailed discussion on its property 
can be found in Ref. 



III. PRETZELOSITY IN THE LIGHT-CONE SU(6) QUARK-DIQUARK MODEL 

The matrix element shown in Eq. flTUl) cannot be solved strictly since a proton state 
contains the nonperturbative information. Conventionally, this state can be expanded in 
a series of light-cone Fock states with coefficients, i.e., the light-cone wave functions, and 
all the wave functions contain the full information of the nucleon. In principle, there is an 
infinite number of wave functions and the wave functions cannot be calculated perturba- 
tively. However, this expansion can be truncated so that only the Fock states with a few 
partons necessary for calculation are left, and the finite number of wave functions can be 
parametrized by experiments or certain models. In this paper, we use the light-cone SU(6) 
quark-diquark model [l^ , in which the proton state is constructed by a valence quark and a 
spectator diquark. The advantage of this model is that some of the gluon and quark effects 
in the spectator can be effectively described by the diquark with a few parameters. However, 
as only valence quark distributions can be directly calculated in this model, one must use 
other information to take into account contributions from sea and gluon distributions. 

Any hadron state can be expanded in terms of a complete set of Fock states at equal 



"light-cone time" [lO|,lll| 

\H) = ^ / [dx][d'^kj_]'(pr,{xi, k±i, \i) 



j-j- u{xiP+, XiP± + fc_Li, A^) -j-j- ejxjP^, XiP± + kjLi, Xj) ^^^^ 
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with the normalization condition 



and the integral over the phase space is 



^ f [dx][d^k^]\M^,,k^,,X,)\'' = 1, (12) 



[dx] =5(1-^2 ^i) n [^^^^] = 167r352(^ k^i) Ylid^k^i/ 1671^). (13) 



i=l 1=1 2=1 i=l 



In the SU(6) quark- diquark model, the proton state with a spin component Sz = ±^ in 
the instant form can be written as 



sin 6(p 
1 



V 



{udfu^ - V2{udyu^ - V2{uu)^d^ + 2{uuYd^] ) + cos9^s{udfu\ (14) 



= — sinOipv {udfu^-V2{ud)-W-V2{uufd^ + 2{uuy^d^ + cos^^5(Mrf) V,(15) 
3 - 

where 6 is the mixing angle that breaks the SU(6) symmetry when 6 ^ tt/A. More explicitly, 
we write the Fock expansion in a diquark model as 

dXqdXD 



^ 167r^ J J 



5(1 -Xq- XD)S^{kg_i + kD±)^p{xj, kjj_, Ai)a]6]|0), 



(16) 



where the wave functions can be extracted from Eq. flMj) . The negative helicity state can 
be written in the same way. 

Notice that now we are working in the instant frame, and we need to transform to 
the light-cone frame. The connection between the two frames is through a Melosh-Wigner 



rotation 



For a spin-i particle, we use qr and qp to denote the instant and light-cone 



spinors, respectively, and the Melosh-Wigner rotation is known to be [1 



T 

If 



UJ 



k^ k~^ + m 




Ml/2 




(17) 



where k^'^ = ki ± ik2, uj = \[xMd + "^q)^ with 



+ 



-, , vTig and 



m£, are mass parameters for the quark and diquark, and Nl^l"^ denotes the Melosh-Wigner 
rotation matrix with the superscript referring to the spin of the particle. For the spin-0 scalar 
diquark, there is no such transformation. For the spin-1 vector diquark, the transformation 
is represented by a 3 x 3 matrix Ai"^, whose explicit expression can also be found in Ref. 18|. 
In practice, we will not use the explicit form, for it does not appear in the final expression 
for the spectator debris due to the unitary property M.^'^M.^ = 1. 

Using Eqs. Q), (ITT]) . (fT6|) . and (ITTI) . we deduce the formula for calculating the trace of 
the correlator: 



E 



1 1 



327r3 xP+ 

j,X,\',XD 

u{xP^, p±, X)Tu{xP+, p±, X') 



^*(x,p_L, A; 1 - X, -p_L, Ad)V^j(x,p_l, A'; 1 



-P±Ad] 



(18) 



TABLE I: Parameters 



niq ms my 
(MeV) (MeV) (MeV) (MeV) 
330 330 600 800 

We substitute io"^+75 for V and notice that |P^) = -^(l-P^) + after a careful calcu- 

lation, we get the result for pretzelosity (the superscript v denotes that it is valid only for 
valence quarks): 

h'[^'''\x,PA_) = X Um'd^lWv. (19) 

oTT'^ y 

with the Melosh-Wigner rotation factor WoiD = V, S) given by 

2M' 

which agrees with the result in Ref. j^, in which a three quark model was used, (pvi^s) 
are the wave functions in the momentum space for the vector(scalar) diquark, which can be 



parametrized by the Brodsky-Huang-Lepage (BHL) prescription [10 



19|: 



M^,p^) = + ^rrf^]}' (21) 



whose normalization should be consistent with Eq. ( 1121) . The parameters such a.s an, the 
quark mass ruq and the diquark mass rriD can be found in Ref. [l^, and we list them in 
Table [H In the same way, the unpolarized distributions can be derived |l5| 

/^^(x,Px) = ^x(lsin'Vv + cos'^V.|), 

/r^(a:,Px) = -i^xlsin'Vv (22) 

Using above formulas, we can calculate pretzelosity and the unpolarized distribution func- 
tion independently, but this model calculation strongly depends on the choice of the wave 
function. The BHL form we use exhibits the sharp fallofF at both the large and small x 
region, but the consequent result for the distribution functions might be inconsistent with 
the realistic situation. For example, the CTEQ6 extraction shows a divergence tendency at 
a; — i> 0, rather than as Eq. (|22|) indicates. 



In order to make our result more close to the realistic situation, we could adopt another 
prescription based on the results of our model. This can be done by combining Eqs. f|T9l) and 
Eq. (I22p . then we get the relation between the unpolarized and pretzelosity distributions: 



Now, we will use the phenomenological extraction of the unpolarized distribution func- 
tions as an input to calculate the pretzelosity distribution. For example, we will use the 
CTEQ6L [2^ parametrization to get fi{x), which have been well tested and constrained 
by many experiments. And then by combining with some assumed transverse momentum 
p_L factor one can get phenomenological fi{x,p±) as input to calculate the pretzelosity 
distribution hj-rp{x,p±). 

Here we need to make some comments on Eq. fl23l) . The left-hand side of the equation 
is a chiral-odd function, which does not have gluon counterparts. But the right-hand side 
is a linear combination of two chiral-even functions, which have gluon counterparts. So the 
two sides have different evolutions, which means that the relation cannot be held exactly at 
any scale. This paradox originates from the fact that the model we use is a no-gluon model, 
which was pointed out in Ref. j^. Nevertheless, we can comprehend this relation from two 
aspects. First, we assume that Eq. (l23l) is valid at an initial scale. Second, the evolution 
effect for pretzelosity appearing on the left is partially (not all) contained in the unpolarized 
distribution. So we conclude that Eq. (!23l) could be approximately satisfied and useful for 
estimating pretzelosity, but we must be careful that the valid scale should not be too large. 

The two approaches to calculate the pretzelosity distribution may lead to different results. 
We denote the former approach, i.e., the model calculation with wave function as input, as 
approach 1, and the latter one, i.e., with phenomenological parametrization of fi{x,p±) as 
input, as approach 2. Here we should emphasize again that we have used a "valence model," 
and we can only directly calculate valence quark distributions in principle. In approach 1, 
for both pretzelosity and unpolarized distribution, only the contribution from valence quarks 
can be considered due to a model calculation. In approach 2, the prescription for pretzelosity 
is the same as that in approach 1, but for unpolarized distribution, the phenomenological 
parametrization can provide us with the information on sea quarks so that we can take 
into account the sea contributions for the unpolarized processes. In this paper, we will use 



ht^-''\x,p^) = [/^^(a:,px) - Ws{x,p^) - ]:fnx,p^)Wy{x,p^), 

2 




(23) 
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both approaches to obtain distribution functions and give predictions separately. Numerical 
results will be presented in Sec. V. 



IV. PRETZELOSITY AND ORBITAL ANGULAR MOMENTUM 



The effect of Melosh-Wigner rotation is also important in other leading twist distribution 
functions such as the helicity and transversity distributions, which have been extensively 
discussed in light-cone formalism in Refs. jisl . 21 1. We find that all the distributions have 
similar form as Eq. (123|) shows, but with different Melosh-Wigner rotation factors. The 



rotation factors for the helicity and transversity distributions are 



{x^+m")^+p'^ —1' respectively. Note the simple relation 



Iq) 



15 



17j and 



D 



so we can immediately yield the equation 



2Ml 



2Ml 



(24) 



(25) 



where gi and hi denote the helicity and transversity distributions, respectively. This relation 
has already been obtained in Ref. [sl with a bag model and a spectator model, and also in 
Ref. jsl with a three-quark model. But in Ref. [3], this relation is not fully satisfied in a 
spectator model where the axial-vector coupling will spoil the relation. We think that this 
issue needs further discussing. 

Now, we should point out that the Melosh-Wigner rotation factor for h^^^'^^ {x,pY) is 
{x^D+mY+p^ ' which is nothing else but the rotation factor Ml(x,p_l) with a minus sign 
introduced in Ref. j^. More importantly, it is the rotation factor for the quark orbital 
angular momentum indicated in Ref. [7], so we see that there is a simple relation between 
the pretzelosity and the quark orbital angular momentum 



-qv. 



-(1)9"/ 



{x,p±) = -h-^T {x,P±) = hl''{x,p±) - gf{x,p±), (26) 

or at the integration level 

L'^^^x) = I d'p^U^x.p^) = -h^i^^'^ix) = hfix) - gT{x). (27) 

We should mention that the orbital angular momentum we denote here is under the definition 
that the quark orbital angular momentum operator is x x (— iV), and the gauge we choose 
is A+ = 0. 
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V. PREDICTIONS ON THE sin(3(/>h -(/>5) ASYMMETRY AT HERMES, COMPASS 
AND JLAB KINEMATICS 



As we pointed out before the pretzelosity has a simple relation with the quark orbital 
angular momentum, thus a measurement of pretzelosity may reveal the information on the 
quark orbital angular momentum. Fortunately, the pretzelosity distribution can be measured 



through sin(30/i — 0s) asymmetry in the SIDIS process 23|, |2J], where the cross section can 
be written as 

^'^^^ - - y + ly')Fuu + sin(30. - 0.)(1 - 



dxdyd(j)sdzd?Phi_ sxy'^ 2 '^'^ 

+other structure functions}, (28) 

with 

Fuu = J'lfiD,], (29) 

FuT =n ^Ml^^ V^i], (30) 

where a compact notation 

T[ufD] = Y.''lj d'p±d'k^S\p^ -k^- Ph^/z) 



ip±,k^)nx,pl)D'^{z,z'kl) (31) 



is used. Then we obtain the sin(30/j — (ps) asymmetry 

2a^ M _ N^sin(3</)h-9is) 
.sin{3(f>h-<f>s) _ sxy^ V » J UT , s 

—,(l-y+-y^)Fuu 

First, we will present the results for the ratio of (the first moment of) pretzelosity (i.e. 
L'^{x) equivalently) and unpolarized distributions, in both approaches. As we mentioned 
above, we can only calculate valence quarks distributions in approach 1, whereas in ap- 
proach 2 we can take into account the sea quark contribution in fi{x). For approach 1, 
we can use Eqs. (fT9l) - (l2T|) directly to calculate. For approach 2, we adopt the CTEQ6L 
parametrization |20| for the unpolarized distribution in this paper and assume a Gaussian 



form factor of transverse momentum as suggested in Ref. 25| : 



A(..Px) = Aw52^^4^. (33) 
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FIG. 1: The ratio hf^^^'^^ (x) / ff (x) . Solid curves for the u quark and dashed curves for the d quark. 
Thin curves correspond to approach 1, while thick curves correspond to approach 2. 



with pI^ = 0.25 GeVl 

Results are shown in Fig.[Tl with a fixed angle 9 = tt/A and an input scale at = 4 GeV^. 
At first sight it seems strange that the ratio of pretzelosity versus unpolarized distribution 
is small in approach 1 compared to that in approach 2, as there is no sea contribution in the 
denominator of approach 1. We can make a brief analysis by taking d quark as an example 
and we could ignore sea quarks for simplicity first. In approach 1, we have 



I 



ff\x) 



l[<fp±^l{x,pl)Wv{x,pi) 
3 



1 |--mj+Pi 



4ai 



+ ^]}Wy{x,pl) 



/rfVexp{-^[^ + ^ 



l-x 



^]} 



- dt e~*Wv{x,4:alx{l - x)t), {t 



4:aj-ix{l — x) ' 



In approach 2, we have 
II = 



ht?'\x) 
ft\x) 



1 / d'^p±dt^''{x,p\)Wv{x,p\) 
3 ^ d^p^dT\x,pl) 

1 
3 



-Wv{x,pI) 



dt e-'Wvix,pij), {t 



p{_. 

Pav 



(34) 



(35) 
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We can easily prove that Wd{x,p\) is an increasing function of and using the inequahty 

4a^x(l -x) <al< pl^, {an = 330MeV, Pav = 500MeV) (36) 

we get that I<II, which means that two ratios can be quite different due to different 
parametrizations. Even if we consider the sea quarks in the denominator, II will be sup- 
pressed, but still can be larger than I. A similar analysis can be applied to the u quark 
though it is a little more complicated, and we can argue that it is possible that II is larger 
than I. From the above analysis we find the different transverse momentum dependence of 
quark distribution functions may cause a big difference in model predictions. This suggests 
the necessity to include transverse momentum dependence in data analysis. Despite the 
differences, both approaches give the prediction that the ratio hf^^^'^^ (x)/ f'Kx) is small, 
which indicates that the relativistic effect, or the quark orbital angular momentum as we 
suggested in our paper, is not significant, if we will integrate the transverse momentum in 
data analysis. 

Before calculating the asymmetry, we still have to know the form for fragmentation 
functions. For the ordinary fragmentation function, we also adopt the Gaussian ansatz 



suggested in Ref. 



25| 



D,iz,z^kl) = D,izf-^^t^^m^ (37) 



TT 



i?2 



with = 0.2 GeV^, and the parametrization for Di{z) can be found in Ref. 26|. The 



Collins function we use in this paper was also given by Anselmino et al. 271], which is also 
a Gaussian function in fact. 

Next, we will present the numerical results of sin(30/j — (ps) asymmetry in the SIDIS 
at different kinematics. For HERMES experiments, only the result on the proton target is 
calculated, while for COMPASS experiments, the proton, deuteron, and neutron targets are 
all assumed. As to Jefferson Lab (JLab) experiments, the result on the proton and neutron^ 
targets are presented. The kinematics are shown in Table [Tll 

The results are shown in Figs. [2]- Fig. HI Obviously, different approaches to distribution 
functions give quite different predictions, and this can be understood from Fig. [H Similarly, 
when we calculate the asymmetry, only valence quarks are summed over in approach 1, 



^ In fact, a '^He target is used to extract neutron data in JLab experiments. 
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TABLE II: Kinematics at HERMES, COMPASS and JLab. 



HERMES 


COMPASS 


JLab 


piab = 27.6 GeV 


piab = 160 GeV 


piab = 12 GeV 


Q2 > 1 GeV^ 


g2 > 1 GeV^ 


g2 > 1 GeV^ 


> 10 GeV^ 


> 25 GeV^ 


W^>A GeV^ 


0.023 < X < 0.4 




0.1 < x < 0.6 


0.1 <y< 0.85 


O.K y < 0.9 


OA<y< 0.85 


0.2 < z < 0.7 


0.2 < z < 1 


OA<z< 0.7 



0.010 

^ 0.005 
w 

-e- 

I 

J 0.000 
*^ -0.005 
-0.010 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 
X 

FIG. 2: sin(3^/i — ^5) asymmetry as a function of x at HERMES kinematics. Solid curves for the 
tt'^ production and dashed curves for the 7r~ production. Thin curves correspond to approach 1, 
while thick curves correspond to approach 2. Only the proton target is assumed here. 

whereas in approach 2, all the flavors are summed over including the sea quarks in the 
denominator. Perhaps approach 2 might be closer to the realistic situation, because the 
parametrization has been well constrained by many experiments and can fit the current 
data much better than a simple model calculation. 

Besides the difference, the two approaches also have something in common. First, we 
find that both approaches predict that the asymmetries are of different sign for tt"*" and 
7r~ productions. Second, we notice that the asymmetry decreases as x increases in the 
valence region, which is quite different from some other known effects, e.g., the sin(0/i + ^5) 
and sm{(f)h — (ps) asymmetries. This can be explained by the ratio of the pretzelosity and 
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FIG. 3: Same as Fig. [2 but at COMPASS kinematics. A for proton target, B for neutron target, 
and C for deuteron target. 

0.010 



0.005 - 



03 

-o- 



J 0.000 
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-0.005 



-0.010 I ^ 



T — r= — ' 1 ' r 




T ' 1 ' 1 ' r 




0.1 0.3 0.5 0.7 0.1 0.3 0.5 0.7 
X X 



FIG. 4: Same as Fig. [2l but at JLab kinematics. A for proton target and B for neutron target. 



unpolarized distribution shown in Fig.dl Another important feature that might be concerned 
by experiments is that the asymmetry is small, up to a maximum less than 1%. Here taking 
approach 2 as an example, we could make a comparison of our result with transversity and 
corresponding sin(0/i + (ps) asymmetry, which have been intensively studied in recent years. 
These two asymmetries can be simply written as 

-(1/2) 

(38) 



, ±(1) rr-L(l/2) 
Asin{3<l>h-<l>s) ^I'l'iT -"l ^sin 



fiD 



The ratio for h^^^ / fi has been shown in Fig. [T], and that for hi/ fi can be found in Ref. [28|. 
We can clearly find that compared with transversity, (the first moment of) pretzelosity is 
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suppressed by 1/2 — 1/10. As to the oscillation function ci;(p_L, A;_l), we could also make an 
estimation on its effect by integrating the angle dependence of the two vectors p± and fc^ 
first, 

r ID 13, . U \ o„2 

(39) 

where we multiply by a weight |-P/i±P just for simplifying the estimation. For a rough 
estimation, we assume {p\) ^ pl^ = 0.25GeV^, (/c^) ^ / {z)"^ = 0.25/(2;)^GeV^, where 
(z) = 0.36 in HERMES. We find that the ratio is about 0.3 ~ 0.4. So from the above 
analysis, the sin(30/i — (ps) asymmetry is suppressed by an order or more compared with 





Id 


(j)pd(j)k\Ph±\^uJi{p±, 




3pl 


C02 


Id 


(j)pd(j)k\Ph±\^(^2{p±, 




2pl + kl 



the sin(0/i + (ps) asymmetry. According to the HERMES data [29|, the size of sm{(f)h + (ps) 
asymmetry is about a few percent, so our result for a small sin{3(f)h — 4>s) asymmetry is 
reasonable from the above expectation. Such small asymmetry seems to be consistent with 



the preliminary COMPASS data 30|, thus it would be a challenge for further experiments 
to measure pretzelosity without transverse momentum information. 

However, the inclusion of transverse momentum dependence in data analysis may provide 
a viable way to extract pretzelosity from sm{3(j)h—4>s) asymmetry measurements. We observe 
that both the Melosh-Wigner rotation factor for hf^^^ and the ratio as Eq. ( l39i) shows are 
increasing functions of p^, which enlightens us that selecting large p± events in data analysis 
might enhance the asymmetry. Unfortunately, the intrinsic parton momentum p_i_ cannot be 
directly manipulated during the measurements. A compromise method is to select large Ph± 
events instead, for Ph± can be directly measured. In this prescription, we can exclude most 
small p± events and ensure that most large p± events will come into data analysis. Next we 
will recalculate our results with a cutoff Ph± > l.OGeV, meanwhile, we will investigate the 
Ph± dependence of the asymmetry. The results are shown in FigJS]- FigE As we expect, 
the X dependence of the asymmetry is indeed enhanced by a few times, thus it might be 
measurable for the experiments, though high statistics would be required due to the cutoff. 
The Ph± dependence of the asymmetry is also presented, and it can reach as large as several 
percent, so it provides an auxiliary measurable quantity to extract pretzelosity. Although 
our proposal to consider to cut off the small Ph± events will enhance the asymmetry, we 
should be aware about this disposal. Our calculation is a leading order approximation and 



based on the factorization for the S 
in the region Aqcd < Ph± < Q 



D 



S, but the factorization was proved to be valid only 



3l| . The ideal kinematic regime to study transverse 
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FIG. 5: The left panel shows the x dependence of the asymmetry with a cutoff Ph± > 1.0 GeV, 
while the right panel shows the Ph± dependence of the asymmetry after integrating all the other 
kinematic variables. The curves used here are the same as that in Fig. [2] 



be important, and a higher order pQCD correction will contribute. 
. 3^, this transition point is around Ph± ~ IGeV, and the parameters 



momentum dependent parton distributions is Ph± ~ Aqcd and not too large Q^. Otherwise, 
the gluon radiation wil 
As pointed out in Ref. 

such as the Gaussian widths we used are suitable at Ph± ^ IGeV. So we must be careful 
about our extension analysis to a larger Ph±, and here we assume that our results at a little 
larger Ph± but not too larger than IGeV are still acceptable. For the experiments, they 
can choose an appropriate cutoff for convenience, and meanwhile we suggest an upper limit 
smaller than 2GeV for Ph±- Another problem is that the events will be strongly suppressed 
at Ph± ^ Aqcd, so it is a challenge for the experiments to collect more data. Anyway, we 
expect further experiments will bring us more information about it. 



VI. SUMMARY 



How to separate the nucleon spin into the intrinsic spin and orbital angular momentum 
parts of the constituents is a fundamental but difficult task. We have shown some debates 
on defining the quark orbital angular momentum in the beginning of our paper and pointed 
out that each definition has its advantage and disadvantage. However, disregarding the 
rationality of the theoretical definition, we can obtain the quark orbital angular momentum 
under each definition, although they might not be the real desired quantity. In Ref. [3], 
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FIG. 6: Same as Fig. O but at COMPASS kinematics. The upper, middle, and lower panels 
correspond to the proton, neutron, and deuteron target, respectively. 

where the expression ip'^x x {—i'V)tp was used as the definition, a quantity L{x) representing 
the quark orbital angular momentum in the light-cone gauge was obtained. In this paper, 
also working in the light-cone representation, we found that L{x) has a simple relation with 
the so-called pretzelosity distribution. 

Pretzelosity is one of the eight leading twist TMD functions, whose feature has been 
discussed in Ref . jsj] . It is known that pretzelosity is considered as the difference between the 
helicity and transversity distributions, reflecting the relativistic effect of the spin structure. 
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FIG. 7: Same as Fig. [H but at JLab kinematics. The upper and lower panels correspond to the 
proton and neutron target, respectively. 

But according to Ref. W, the moment of pretzelosity is nothing but the quark orbital 
angular momentum. So it provides a new possibility to access the quark orbital angular 
momentum through pretzelosity, and fortunately, pretzelosity can be measured in SIDIS 
through sin(30/i — (ps) asymmetry. In our paper, under the framework of the SU(6) quark- 
diquark model, we use two approaches to calculate pretzelosity, and we present our numerical 
prediction on the sm{3(f)h—4's) asymmetry under three kinematics for HERMES, COMPASS, 
and JLab experiments. If we do not apply any extra constrained condition, our prediction 
shows a small asymmetry that might bring a great challenge for experiments. However, 
after applying a cutoff on Ph±, we get an enhanced asymmetry which is measurable for 
experiments, though we should still be careful that Ph± must not be too large. We expect 
that future measurements can bring us exciting new insights on the spin of the nucleon by 
including information of transverse momentum dependence in data analysis. 
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